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Abstract. We compare metrics obtained through analytic perturbation theory with their numerical 
counterparts. The analytic solutions are computed with the CMMR post-Minkowskian and slow 
rotation approximation due to Cabezas et al. [1] for an asymptotically flat stationary spacetime 
containing a rotating perfect fluid compact source. The same spacetime is studied with the AKM 
numerical multi-domain spectral code [2, 3] . We then study their differences inside the source, 
near the infinity and in the matching surface, or equivalently, the global character of the analytic 
On perturbation scheme. 

CN 

Keywords: Relativistic astrophysics, post-Minkowskian approximation. Harmonic coordinates, 
I— I Rotating stars 

O PACS: 04.25.Nx, 04.40.Dg 

^ INTRODUCTION 

^ Despite the great effort invested, there is still no exact solution of Einstein's equations 

\0 able to describe a stellar model, i.e., a singularity-free rotating body that has been 

matched to an asymptotically flat vacuum exterior. In the last two decades, the attention 
has moved to the field of approximate solutions. Among the recent ones is the AKM 
^ code [2, 3]. It is a multi-domain spectral method, and the difficulties many other codes 

O have on surfaces of discontinuity of some sources due to Gibbs phenomena are solved 

computing the solution of the different domains and then imposing matching conditions. 
T1 The number n of Chebyshev polynomials in the expansions, resolution and equation of 

. ^ state (EOS) can be chosen, and reaches machine accuracy for high enough n. We will 

^ use it to check the behaviour of the CMMR post-Minkowskian+slow rotation analytic 

^ approximation scheme some of us introduced in [1]. We begin fixing the properties of 

the problem spacetime while briefing the CMMR basics. Then we study the relative error 
in the metric functions between schemes and some physical properties of the source. 



BUILDING THE METRICS 

The spacetime studied Y is stationary, with timelike Killing vector field ^ , and axisym- 
metric, being t] the associated closed-orbits spacelike Killing vector field that satisfies 
regularity condition near the axis. It is built from the matching of two spacetimes. The 
first one, is filled with a perfect fluid in circular flow so that its velocity can be writ- 
ten u = + (or\), with 1// adjusted to make u^Ua = — 1. The function (O is constant, 
making the fluid rigidly rotating. It has constant energy density, /i = /io, so integrating 
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FIGURE 1. Relative error between CMMR and AKM in (a-b) gtt and gtcp, Mb = ^ x 10"^ CO = 0.2 
(A ^ 0.0030, a ^ 0.098); (c-d) and gtcp, Mb = ^x 10-^, co = 0.2 (A 0.013, Q. ^ 0.098). Note that 
both the scale and the angular dependence of the error decrease with bigger Mb. The thin dotted lines 
represent the AKM and CMMR surfaces (indistinguishable in this picture size). 



Euler's equations gives the pressure p — fio [(i//^/ i/^e) — 1], with xj/ = xj/^ for p = 0. The 
second spacetime, Y^, is asymptotically flat vacuum surrounding y~. It is not restric- 
tive for us to use global harmonic Cartesian-like coordinates {x^} in y[4]. Working 
in spherical-like coordinates associated to {x^}, the = surface E on which y~ and 
are matched can be written as an expansion r^: = [l + CJ Q?P2{cos 6)] + ^{Q.^) in 
Legendre polynomials P^, with a a constant. It has been truncated introducing a slow ro- 
tation approximation parameter Q we have chosen as = (O^r^/m, where m = ^Tiiior^ 
is the Newtonian mass of a sphere of radius r^. 

To solve Einstein's equations, we use a multipolar post-Minkowskian approximation 
as follows. Defining a parameter A = m/r^, the exact metric in each spacetime 
is decomposed as g^(A,^2) = T] + h=^(A,^2), with rj the flat metric. Then, Einstein's 
equations are solved iteratively in A. Both h^(A,^2) are tensor spherical harmonic 
expansions that are truncated, in this case, to contain powers lower than Q!^. This 
restricts the number of multipole moments M/, appearing in the exterior solution. 
We then match y~ and imposing continuity of the metric and its first derivatives. 
This fixes all coefficients and the stellar model depends then only on /io, CO and r^. 



The AKM code computes the matched value of the functions U,k,W and a in the 
general line element of a stationary axisymmetric perfect fluid or asymptotically flat 
vacuum spacetime in quasi-isotropic coordinates (see, e.g. [2], where {p, ^} are cylin- 
drical associated to quasi-isotropic coordinates {r, G}) at each point of a coordinate grid 
of user-definable resolution. It also gives a lot of information in terms of physical and 
geometric parameters, such as multipolar moments Mq and J\, baryonic mass Mb, an- 
gular velocity 0), equatorial radius and central pressure pc among others. Once the 
values of two of them and the EOS have been fixed, the code can compute the metric. 

COMPARISON RESULTS 

To compare the results of CMMR and AKM for /i = /io, we must first find the change 
of coordinates from the spherical-like ones of CMMR to quasi-isotropic. This change 
is necessarily approximate, introducing a new source of error in the comparison. This 
makes the relative error we compute between the metric functions of each scheme at 
a point to be a strict upper bound. For the comparison, CMMR was computed up 
to order {^{?i^/^),^{Q?)), and AKM was set to use 12 Chebyshev polynomials in 
each direction. Then, working in dimensionless quantities (G = c = ^Uq = 1) we must 
choose which two parameters to fix in both CMMR and AKM. For this work, we have 
dealt with two sets, first {Mq, co} and then {r^, co}. Once co is fixed, CMMR results 
depend only on r^. We get its value equating both Mq (alternatively, r^) values. The 
Mo adjustment gives better results and is the one we will focus on. Figs. 1-2 show the 
relative error in gtt and gt(p (gu plots are very similar to gtt ones) on a quadrant of the 
plane p - C f or O) = 0.2 and AKM values of = 8 x 10"^ 8 x IQ-^ and 8 x 10"^ 
For a typical neutron star density /io = 4 x lO^^kgm"^, they would correspond to a 
frequency v ^ 1033 s"^ and Mo ^ O.OO3M0, O.O3M0 and O.3M0, respectively. Table 1 
shows their CMMR values and relative errors of some quantities. The rather extreme 
cases of 0) = 0.7 are included to check the behaviour of our slow rotation approximation 




0.00 0.05 0.10 0.15 0.20 0.00 0.05 0.10 0.15 0.20 

P P 

FIGURE 2. Relative error between CMMR and AKM in gtt and gt(p for Mb = S x 10"^ (O = 0.2 
(A ^ 0.056, a ^ 0.098). (Continued from Fig. 1) 



TABLE 1. CMMR values of some quantities and relative error with AKM £ = 

for two members ((0 = 0.2 and CO = 0.7) from each studied sequence: Mb = 8 x 10~^, 8 x 10~^ 
and 8 x 10"^ 





CMMR 




CMMR 


£ 


CMMR 


£ 


CO 


0.2 




0.7 




0.2 




Mo 


0.00007985 




0.00007986 




0.00079334 




Jl 


4.6065620 e-9 


0.00020 


1.75490113 e-8 


0.017 


2.13041342 e-7 


0.0021 


re 


0.02674082 


4.4 e-6 


0.02792957 


0.0074 


0.05689761 


0.00073 


Pc 


0.00147636 


0.012 


0.00136990 


0.00041 


0.00668441 


0.056 


rs 


0.02663511 




0.02663467 




0.05667470 




CO 


0.7 




0.2 




0.7 




Mo 


0.00079362 




0.00769003 




0.00770289 




Jl 


8.10712355 e-7 


0.018 


9.52547199 e-6 


0.032 


0.00003614 


0.043 


re 


0.05940141 


0.0065 


0.11664607 


0.014 


0.12158066 


0.0080 


Pc 


0.00620181 


0.044 


0.02810237 


0.27 


0.02607038 


0.26 


rs 


0.05667105 




0.11620612 




0.11619185 





for high values of ^2 (^2 ^ 0.49ft)). 

For the three cases studied, relative errors in metric functions increase roughly two 
orders of magnitude if we make Mb ten times bigger, being higher for gt(p and the 
interiors. This is expectable since we have fixed Mq, i.e. the behaviour of gn near spatial 
infinity, what can cause the high values of £{pc) as well. Significant error discontinuities 
are located on equatorial/polar lobes (Fig. lb) and can be caused by the truncation at 
&{Q?) of and the Legendre expansion of £. This is supported by the smooth plots 
we get in the static limit with adjustment, and the increased lobular appearance when 
Mb decreases (giving rise to more oblate configurations for the same ft)). We expect this 
angular dependence of e to decrease including more terms of the series. 

The error inside the source is systematically bigger than outside it, but comparable. 
We plan to use pc and J\ to fix and expect better results in the interior. We will add 
new physical quantities to the comparison to see how much the general performance 
improves going further in the approximation as well as other EOS. 
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